Strong lensing of gravitational waves as seen by LISA 
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We discuss strong gravitational lensing of gravitational waves from merging of massive black hole binaries 
in the context of the LISA mission. Detection of multiple events would provide invaluable information on 
competing theories of gravity, evolution and formation of structures and, with complementary observations, 
constraints on Hq and other cosmological parameters. Most of the optical depth for lensing is provided by 
intervening massive galactic halos, for which wave optics effects are negligible. Probabilities to observe multiple 
events are sizable for a broad range of formation histories. For the most optimistic models, up to < 4 multiple 
events with a signal to noise ratio > 8 are expected in a 5-year mission. Chances are significant even for 
conservative models with either light (< 60%) or heavy (< 40%) seeds. Due to lensing amplification, some 
intrinsically too faint signals are brought over threshold (< 2 per year). 

PACS numbers: 95.30.Sf, 95.85.Sz 
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The space-based Laser Interferometer Space Antenna yj, 
LISA] is expected to observe gravitational waves (GWs) and 
open a new window for astronomy. Coalescencing massive 
black hole binaries (MBHBs) with total masses in the range 
10 3 -10 7 Mq out to z ~ 10-15 are expected to provide the 
loudest GW signals at LISA frequencies, / ~ mHz In 
the standard cold dark matter (CDM) hierarchical cosmology, 
MBHBs form in large number during the multiple mergers 
experienced by their host galaxies. LISA should detect from 
few to several hundreds coalescences per year. Within such 
promising premises we discuss a potential new chapter for 
LISA science: multiple imaging of distant sources by inter- 
vening lensing galaxies. Gravitational lensing statistics, in ei- 
ther quasar or radio-galaxy surveys, is usually considered on 
a sample of several thousands of sources. Due to the unprece- 
dented high redshift of LISA sources, and the related very 
high optical depth for lensing, multiple events are possible 
even for hundreds of detections. Astrophysical and theoreti- 
cal rewards might be valuable, i) Lensing of GWs would be 
another impressive confirmation of the theory of general rel- 
ativity. How GWs propagate near a massive body might shed 
light on competing theories of gravity, ii) Constraints on cos- 
mological parameters might be obtained in the range z > 10. 
Measurements of time-delay, which can be accurately deter- 
mined for transients [3], might be very useful. Hi) Lensing 
statistics might inform of the growth and structure of mass 
halos at z < 3. The expected number of events strongly de- 
pends on the form and evolution of the galaxy number density. 
iv) The magnification effect could help in finding electromag- 
netic counterparts and observing objects otherwise too distant 
or too faint. 

Strong lensing by ground-based GW detectors was dis- 
cussed in J4(]. Here, we consider a lower frequency range 
where wave optics could play a role. Diffraction in lensing 
of GW is well understood [5] and is effective only for small 



lenses [6]. Our study on strong lensing is complementary to 
those on the weak lensing distortion of GW signals, which 
mainly focused on their use as standard sirens 0]. By default, 
we assume a flat ACDM model with Qm = 1 — Q\ = 0.3, 
and H = 70 km s _1 Mpc _1 . 

Lenses. Statistics of strong lenses is a well assessed astro- 
nomical tool |3|. Here, we mainly follow lalSI]- The differen- 
tial probability of a source to be lensed is 
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where a and Zd are the velocity dispersion and the redshift of 
the deflector, respectively, s cr is the cross section and dn/ da 
is the differential lens number density, dn/da can be modeled 
as a modified Schechter function 11011 . This is accurate up to 
Zd < 1, where galaxies provide the bulk of the cross section 
for LISA sources. We took a conservative approach. Together 
with a constant comoving number density, we also considered 
a pessimistic scenario, i.e., an evolving case with a smaller 
number of lenses at high z 111 111 . 

We account only for lensing by early type galaxies, mod- 
elled as singular isothermal spheres (SISs). Their Einstein 
radius is Re = 47r(cr / c) 2 DdDds/ D s ; Dd s , Aj and D s are 
the angular diameter distances between the deflector and the 
source, and the observer and the lens or the source, respec- 
tively. Two images form at x± = y ± 1 if y < 1, with flux 
magnification /u± = (1/y) ± 1; x and y are the image and 
the source position normalized to the Einstein angular radius. 
The GW form is amplified by A± = */ju±. y can be mea- 
sured through the ratio of the A±. The delay between the 
arrival time of the images, At = t_ — t+, is 



At = At z y, At z = 
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For <7 ~ 200 km s , Zd — 5 and z s ci 
LISA angular resolution is quite poor, 
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FIG. 1. y max as a function of for z s = 10. Thick and thin lines 
are for t+ — 2 or 4 years from the survey kick-off. 



expected to be very massive and luminous. In case the deflec- 
tor is located, one can precisely constrain the position of the 
source and try to detect it with follow-up observations. Com- 
bining measurements of At, z<j and z s , and D a from the GW 
analysis might allow us to constrain Ho and dark energy. A 
cosmographic approach independent of time-delay is based on 
lensing statistics 13J. 

The cross section for lensing is s cr = ^R^iVmax — 2/min)> 
where y m i„ < y < y ma x is the allowed range to form de- 
tectable multiple images. y max depends on the lens mass and 
redshift, the threshold signal to noise ratio (SNR t h), the un- 
lensed amplitude (SNRj nt ), see Eq. [3] the arrival time (t+), 
and the total survey time (T sur , see Fig. [TJ. We considered 
?sur = 5 years. The emission frequency enters in y m i n , which 
excises the region near the central caustic where wave optics 
is effective. Geometric optics is valid for y > y m i n . 

Lens discovery rates are affected by the ability to resolve 
multiple images [3]. Amplification has to push the signal 



above threshold, A± 
source position to 



> SNR^/SNRi u t, which limits the 



y± < y„ 
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according as the requirement is on the + or the — image, re- 
spectively. Multiple events are detectable if y < y^ ax , which 
goes to 1 only for SNR int > SNR t ~ . 

Due to the finite duration of the survey, statistics of transient 
phenomena involve some missing events due to time delay 
03J. To observe both images, we require that At < T SUT — t + , 
which further constrains the source position. In Fig. Q] we 
plot y max due to the finite observation time as a function of 
z&. This constraint plays a role only for very massive lenses, 
which are very rare, or for very late arrival times. 

Wave-optics. Since we cannot observe multiple images 
when the interference pattern is pronounced, we have to dis- 
card configurations where wave effects are large Q3Q . Diffrac- 
tion is negligible when the wavelength is much smaller than 
the gravitational radius of the lens, i.e., At ^S> f~ l ta|. 
Geometric optics is valid if the coherence time of the sig- 
nal is much smaller than the time delay between images, 
A/At ^> 1, where A/ is either the bandwidth of the detector 



FIG. 2. Relative error in the GW amplitude introduced by ne- 
glecting diffraction. Contour values are 0.2, 0.1 and 0.05; darker 
areas correspond to larger errors. Time is in units of At z . Above the 
full line, At x / > 10, geometric optics is valid. 
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FIG. 3. Contour plot of y m i n (such that At x / > 10) for z<j = 2 and 
z s = 6. Full lines take values from 0.1 to 1 in steps of 0.1; dashed 
lines are for y m i n = 10 -2 , 10 -3 and 10 -4 . 



or the range in emission frequency. These conditions deter- 
mine y m i n . The error in the estimated wave amplitude intro- 
duced by neglecting diffraction is displayed in Fig. |2j where 
we exploited the sum representation of the diffraction integral 
fl. 

Two features determine the effective frequency for MB- 
HBs. First, most of the signal is emitted near the coales- 
cence time. Second, LISA sensitivity drops for / > / c h = 
3 x 10~ 3 Hz. The effective frequency is then the minimum 
between / at coalescence and / c h- Since the final frequency 
is much larger than the initial one, the two conditions give 
very similar constraints. Apart from a very small region near 
the central caustic, see Fig. [3] wave effects are negligible. 
For a source z s = 6 emitting at / c h behind a massive lens 
(a ~ 150 km/s, z d = 2), y rain ~ lO" 2 . 

Optical depth. The lensing probability for a source is the to- 
tal optical depth r, obtained by integrating Eq. (HJ; it is plotted 
in Fig.|4]as a function of z$. In our standard case, J7m = 0.3, 
the comoving dn/da is constant and SNR^ = 8 is the de- 
tection threshold. The more effective constraint on the cross 
section is the requirement of a loud enough second image. If 
we either lower SNR^ or consider intrinsically brighter sig- 
nals, r dramatically increases, r is very sensitive to the value 
of A and goes down for a smaller number of lenses at high z. 



3 




FIG. 6. Probability of detecting N lensing events in the standard 
Z s case; thick curves include 68% confidence regions; "L" = Light. 



FIG. 4. t for / = / c h and t+ = 2 years. Thin and thick lines are for 
SNRi n t — 6 or 20, respectively. Full lines are for the standard case. 




FIG. 5. r contour plot for a MBHM with / = / c h and t+ — 2 years. 
Thin and thick lines are for M2/M1 = 0.1 or 1, respectively. Values 
go from to 4.5 x 10~ 3 in steps of 0.5 x 10~ 3 . 



Probabilities are quite large for a range of BH binaries with 
masses (10 4 -10 7 M Q ) and redshifts (z < 15) of astrophysical 
relevance in the LISA window, see Fig. [5] r is plotted in the 
standard case as a function of z and the main BH mass, Mi, 
for two different mass ratios. The SNRj nt of a given merger 
was computed according to [2]. Given M\, the larger the mass 
ratio, the larger r. 

LISA sources. To estimate the lensing probabilities we 
need plausible merger rates. We adopt MBH assembly his- 
tories created via dedicated Monte Carlo merger tree simula- 
tions framed in the hierarchical paradigm 01211 . Either first 
seed BHs are light, M > IOOMq, being the remnant of 



the first POPIII star explosions il 1 3fl or already quite heavy 
(M > 10 4 M Q ) seeds form by direct collapse of proto-galactic 
discs II 1 411 . or through efficient accretion onto a "quasi star" 
forming in dense disks prone to bars-within-bars instabilities 
1 15]. The first two models we analyzed (Light and Heavy 
models, LE and SE in I1 1611 . respectively) have been adopted 



as a benchmark for model comparison by the LISA parame- 
ter estimation taskforce. Furthermore, we considered a pair of 
"hybrid" models where the two formation processes coexist 
and the seed mass spectrum spans the range 10 2 -10 5 Mq lfl7ll . 
The two hybrid models have somewhat different formation ef- 
ficiencies, which are set by the maximum spin parameter Axhr 
of the halo that allows an efficient inflow of gas in the nuclear 
region. We considered a conservative version (Atiu = 0.01, 
Hybrid!) and a more optimistic one (Axhr = 0.02, Hybrid!!). 

All models have been extensively tested against observa- 
tional constraints, such as the present day mass density of nu- 
clear MBHs, the optical and X-ray luminosity functions of 
quasars and the unresolved X-ray background 11211 . Light, 
Heavy and Hybrid! models are on the conservative side, while 
Hybrid!! is more optimistic. LISA source redshifts are much 
higher than usual, see TableU and allow to investigate the Uni- 
verse to very large distances. By comparison, radio-sources in 
the CLASS survey had z s = 1.4 ± 0.7 0]. 

Expected events. The mean number of lenses expected to be 
observed, (N), is given by the sum of the optical depths of all 
the sources. We simulated the intrinsic properties instead of 
the observed distribution of amplitudes so that we did not have 
to correct for any magnification bias yfl. Results are listed 
in Table U The Poissonian probability for multiple events is 
plotted in Fig. [6] Chances are fairly large for all the build-up 
scenarios. Wave effects are negligible: in the geometric optics 
limit (y m in = 0), r tot increases by < 0.01%. 

The bulk of the total cross section is due to high SNRj nt 
events. In the Light scenario, the fraction of events (~ 54%) 
intrinsically below threshold contributes only < 8% of Tt t- 
Due to the many events with very high SNRj nt , results for 
the Heavy case are not affected at all by lowering SNR^, 
whereas probabilities increase significantly in the Light and 
Hybrid! cases. LISA should detect independently two above 
threshold signals in the same position in order to claim lens- 
ing. Any dedicated pipeline for detection of double peaked 
events might strong enhance lensing detection. 

A larger, but still compatible with observations, A strongly 
increases the probabilities. New merger trees should be gen- 
erated, but to stress the role played by a larger A we focused 
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scenario Light (90, z s = 10.4 ± 4.1) Heavy (30, z s = 5.8 ± 2.4) Hybridl (130, z s = 7.2 ± 3.8) Hybridll (390, z s = 6.0 ± 3.4) 



& SNR" h 


{N) N > 1 N > 2 (N) N > 1 N > 2 (JV) AT > 1 N > 2 (N) N > 1 AT > 2 


0.3 co. 8 


.38 ±.03 32 ±2 5. 8 ±.7 


.24 ±.02 21 ±1 2.4 ±.3 


.55 ±.04 42 ±2 10.6 ±1.2 


2.83 ±.09 94.1 ±.5 77.4 ±1.5 


0.2 co. 8 


.58±.04 44 ±2 11 ±1 


.35 ±.03 29 ±2 4.8 ±.6 


.81 ±.06 55 ±3 19 ±2 


4.17 ±.14 98.4 ±.2 92.0 ±.9 


0.3 ev. 8 


.30 ±.02 26 ± 1 3.7 ±.4 


.19 ± .01 17±1 1.5 ±.2 


.44 ±.03 35 ±2 7.2 ±.9 


2.24 ±.07 89.4 ±.7 65.7 ± 1.7 


0.3 co. 5 


.52 ±.03 41 ±2 10 ± 1 


.24 ±.02 21 ±1 2.5 ±.3 


.68 ±.05 49 ±2 14.7 ±1.6 


3.06 ±.10 95.3 ±.5 80.9 ± 1.4 



TABLE I. Lensing probabilities assuming SNR^ > 8. (TV) and chances (in %) to observe one or more lensing events, TV > 1, 2, for different 
hypotheses (name, number of events per year, typical source redshift and its dispersion). Biweight estimators are reported. 



on the variation of r. Reducing the number of lenses at high 
z has an opposite effect but chances are still sizable. 

Lensing amplification can make some intrinsically too faint 
mergers luminous enough to be detected. If we are not inter- 
ested in the second image, y max can be larger than one and is 
determined only by the flux condition on the ± image. In the 
Light {Hybridl) model, nearly 3(2)% of the otherwise unde- 
tected events, > 1 event per year, may become loud enough to 
be seen. In the Hybridll scenario, 8 ± 3 additional events per 
5 years are expected. 

We provided one of the first attempts at formalizing lensing 
statistics in presence of wave-effects. Our assumptions were 
conservative: a 5 years mission; high detection thresholds; 
not so optimistic formation scenarios; no account of lensing 
by late-type galaxies, whose distribution is poorly known but 
contribute < 30% of T to t A- Probabilities to observe multi- 
ple events were anyway sizable, from < 20 to < 100%. Lens- 
ing statistics depend on the cosmological parameters, suggest- 
ing new potential tools for cosmography with LISA. In case 
of optical identification of lenses, additional measurements of 
time-delay might help to constrain Hq and dark energy, r is 
very sensitive to dn/da too, which might help to understand 
evolution of early-type galaxies up to high redshifts. 

The main unsureness in our forecast is due to the uncertain 
build-up process. Tight predictions need knowledge of galaxy 
evolution at high z, one of the main LISA goals. To establish 
that lensing probabilities for LISA might be significant and to 
point to the relative rewards, it was enough to consider a broad 
range of formation scenarios without trying to maximize Ttot- 
Hopefully, LISA lensing might help to constrain the formation 
history. 

Apart from astronomical applications, the detection of even 
a single GW lensing event would provide unique information 
on competing theories of gravity by allowing to test velocity 
and propagation of gravity. Up to now, GWs are still missing 
on an experimental ground. Their mere detection by LISA 
would be a success. Serendipitous GW lensing might be an 
extremely beneficial bonus. 
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